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Abstract. Binary merging constitutes a complementary 
mode of galaxy evolution to the canonical hierarchical 
clustering theory. This merger-driven evolution not only 
influences the galaxy mass distribution function, but may 
drive the main galactic activity cycles, such as star bursts 
and the activity in the nuclei (AGNs). In this paper, we 
use galaxy aggregation dynamics together with a possible 
merger-driven starburst and AGN phenomena, to study 
the effects of evolution of these consequential activities 
for the faint IRAS source counts and the infra-red back- 
ground. We find that the strong evolution of IRAS 60 fim 
source count at flux range of 10 m Jy ~ 1 Jy is difficult to 
be explained only by the merger rate decrease with cosmic 
time. We need to assume a redshift dependent infrared 
burst phase of ultraluminous infrared galaxies (ULIGs) 
from the gas rich mergers at high redshift and the gas 
poor mergers at low redshift. The background intensity at 
60 jim which we get from this aggregation evolution is a 



lower limit of only 1.9 nWm 



' 1 , about half of those 



estimated by some previous models, but close to the lower 
end of the range derived by Malkan & Stecker (1998). 

Key words: Galaxies: starburst - Galaxies: Seyfert - 
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actions - Galaxies: luminosity function, mass function 



1. Introduction 

Over the last few years, observational progress in discov- 
ering the faint Universe at ultraviolet (UV)/optical and 
infrared wavelength via the groundbased redshift surveys, 
Hubble Space Telescope (HST) multi-color images and the 
infrared satellites such as IRAS (Infrared Astronomical 
Satellite), ISO (Infrared Space Observatory) and COBE 
(Cosmic Background Explorer), challenges and renews our 
present understanding of cosmological structure formation 
and the evolution of luminous galaxies, which are usually 
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based on two competing scenarios "traditional monolithic 
collapsing" versus "hierarchical clustering theories". The 
main difference between these two kinds of models is the 
formation mechanism and the first epoch of the formation 
of massive galaxies. 

More recently it has been shown that the emission his- 
tory (e.g. Madau plot) derived from UV/optical detections 
agrees quite well with the star formation history indicated 
by hierarchical clustering theories (Madau 1997, Madau et 
al. 1996, 1998, Ellis et al. 1996). Since the FIR/submm un- 
veiled a dust-shrouded early active phase and the recent 
detection of high redshift evolving galaxies by deep NIC- 
MOS/VLT images (Lilly et al. 1998, Benitez et al. 1998), 
the discussion of the structure formation and the scheme 
of galaxy evolution is going to require more subtlety and 
sophistication. For a clear picture of the star formation 
history and the galaxy evolution, we probably need more 
information at high redshift, especially from FIR/submm 
deep surveys. 

A successful interpretation of the excess of faint blue 
galaxies within the merger-trigger starburst scenario in 
a hierarchical clustering universe has been presented by 
Cavaliere & Menci (1997) using binary aggregation dy- 
namics. 

The binary aggregation dynamical approach for the 
galaxy evolution presented by Cavaliere & Menci (1993, 
1997), which includes more dynamics to describe a fur- 
ther step in galaxy-galaxy interactions within the scheme 
of direct hierarchical clustering (DHCs) , probably can help 
to alleviate some intrinsic problems in the DHC scenario, 
such as the overproduction of small objects and the dif- 
ficulty of a reconciliation between the excess faint counts 
and the flat local luminosity function. While, on the other 
hand, binary merging which plays a different and com- 
plementary role in structure formation, can actually con- 
tinue to flatten the mass distribution function N(m, z) 
and hence the luminosity function N(L, z) till even mod- 
erate redshift z < 1. 

Although there are many other possible evolutionary 
scenarios which could interpret the present observations 
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in this or that way, the reason that we are encouraged 
to explore here a merger-driven galaxy evolution picture 
with the binary aggregation dynamics, is simply because 
the IRAS database does show that most of the luminous 
infrared sources are actually interacting/merging systems 
(Kleinmann & Keel 1987, Sanders & Mirabel 1996, Hutch- 
ings & Neff 1987, Vader & Simon 1987a,b), and many 
merger events can even occur till moderate redshift (z < 

I)- 

In this paper, starting from the basic picture that 
galaxy evolution is driven by the galaxy-galaxy interac- 
tion complementary to the standard hierarchical cluster- 
ing scenario (DHCs), we investigate a simple model with 
the basic concepts: 1) galaxy-galaxy interaction can even 
occur till quite moderate redshifts [z < 1) within large 
scale structures that could offer the best combination of 
volume and density contrast, and it could erase the dwarf 
part of the mass distribution and produce the massive tail, 
thus flatten the luminosity functions; 2) starburst/ AGN 
activities may be triggered by the merger events during 
the structure formation and evolution. This is actually 
the point different from other models where only the star- 
burst is considered as the consequence of merger events 
(Cavaliere & Menci 1997, Somerville & Primack 1998, 
Somerville et al. 1998); 3) the gas rich merger events at 
high redshift may trigger drastic starburst and AGN activ- 
ities in the central region of the merging galaxies and thus 
enhance dramatically the FIR luminosity because of the 
accumulating of dust grains from their progenitor galax- 
ies as well as the dust newly formed in the star forming 
regions. This FIR luminosity burst phase is usually called 
ultraluminous infrared galaxies (ULIGs); 4) in our model, 
we assume a redshift dependent infrared burst phase for 
the gas rich merger events at high redshift and the gas 
poor mergers at low redshift, which means the enhance- 
ment of the infrared luminosities by the high redshift gas 
rich mergers is much higher than those of low redshift gas 
poor mergers. The motivation for this kind of thinking 
is: 1) direct observation of enhanced starburst-AGN ac- 
tivities in interacting galaxies, especially the ULIGs. This 
extremely infrared bright burst phase is believed due to 
the starburst merger events with far infrared luminosity 
L fi r enhanced both by the accumulation of the dust mass 
Md and the increase of dust temperation Td as the rela- 
tion Lfi r oc Md T|. This burst phase could increase the 
infrared luminosity by a factor of about 20 over that of 
normal starburst galaxies. (Ashby et al. 1992, Terlevich 
& Melnick 1985, Heckman 1997, Perry & Williams 1993, 
Taniguchi & Ohyama 1998); 2) numerical simulation of 
the starburst/ AGN evolutions by galaxy mergers, espe- 
cially for their correlation and the burst phase (Wang 
& Biermann 1998, 1999, Wang 1999); 3) the successful 
interpretation of the excess faint blue counts by galaxy 
merging and the consequential starburst activities (Rocca- 
Volmerange & Guiderdoni 1990, Cavaliere & Menci 1997, 
Carlberg 1992, Calberg k Chariot 1992); 4) the pair pro- 



duction absorption of high energy 7 rays by intergalactic 
low energy photons is expected to produce a high energy 
cutoff in extragalactic 7 ray spectra. The new data of 
Mkn501 appear to show an extension of the TeV 7 ray 
spectrum till about 24 TeV, which sets an upper limit 
for the intergalactic infrared emission history during the 
structure formation (Aharonian et al. 1999). 

Since the dust-shrouded geometry could strongly af- 
fect the AGN spectrum and cause significant radiation 
at infrared wavelength in the framework of a Unification 
Scheme, we thus add the statistics of nuclear activity in 
galaxies as one additional constraint in the context of the 
various models already in the literature. We know that nu- 
clear activity is detected to very high redshift and provides 
a strong constraint on the cosmological evolution models. 
Therefore we are less conservative than Malkan & Stecker 
(1998), but use more constraints than, e.g., Somerville & 
Primack (1998). Obviously, any such modelling depends 
on the basic concepts used, and thus our results provide a 
strong lower limit to the far-infrared background. 

The outline of this paper is as follows: 1) in Sect. |[ 
we introduce the binary aggregation dynamics, and the 
models of the interaction kernel; 2) we discuss the pre- 
scriptions of mass-to-light ratio in our model for the lu- 
minous infrared galaxies which is enhanced by the star- 
burst and AGN activities in Sect. [| especially modelling 
of the redshift dependent infrared burst phase for ultra- 
luminous infrared galaxies (ULIGs) from gas rich mergers 
at high redshift and a suppressed burst phase for gas poor 
mergers at low redshift; 3) in Sect. ||, we discuss the nu- 
merical simulations and compare the Monte Carlo results 
with the IRAS 60 fim source counts from three major in- 
frared sources (starburst galaxies, Seyferts and spirals). 
We check also the redshift distribution of these infrared 
sources which are brighter than Sqq ~ lOmJy, calculate 
the integrated background level at wavelength of 60 /im, 
and make an extrapolation to 25 [im, and 100 fj,m based 
on the model spectrum and souce counts. Finally, we give 
our conclusions. 

2. Binary aggregation theory 

The classical approach to aggregation phenomena is based 
on the Smoluchowski(1916) equation. We start with the 
continuous form: 



<9N(m,t) 
di 

m— m* 



dm K(m , m — m . t)N(m , t)N(m — m , t) 



-N(m,t) I dm' K(m,m ,t)N(m ,t) (1) 



where N(m,t) is the mass distribution function in "co- 
moving" form, which describes the number density of 
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galaxies within the mass range (to, m + dm) at cosmic 
time t . The mass to = MjM* is normalized in terms of M* , 
which is the mass corresponding locally to the standard 
characteristic luminosity L* (see Peebles 1993). to* and 
to** represent the lower and upper limits on the masses 
of galaxies. Usually we can set to* = 0, to** = oo in case 
when to* is actually much smaller than the total mass of 
considered system. 

The kernel K(m, m ,t) = n g {t) < > reflects the 

interaction rate for each pair of masses (to, to ), and de- 
pends on the mechanism and environment of such aggre- 
gations. n g is the density of galaxies, it scales on average 
with the expansion of Universe as n g cx (1 + z) 3 ; V(t) is 
the relative velocity of the interacting pairs, and £ is the 
cross section. The average of them runs over the galaxy 
velocity distribution. 

In aggregation dynamics, the interaction kernel K = 
n g < TjV > is the key point for driving the whole evolu- 
tionary process, which depends strongly on the environ- 
ment structures. 

In this context, we discuss simply the case when the 
colliding galaxies are in certain structures, thus the cross 
section could be assumed to be the hyperbolic encounter 
prescription by Saslaw (1985), Cavaliere et al. (1991, 
1992), which is in the form: 



G(m + to ) 
r m V 2 



(2) 



with r m — r^m 1 / 3 ,v m ~ v*m 1//3 , where t>* and r* are the 
velocity dispersion and the dark halo radius of the present 
M* galaxy. The function e ^~^j describes the decreasing 
efficiency of the aggregations with an increasing relative 
velocities, which could be determined by the N-body sim- 
ulations (see, Richstone & Malumuth 1983). 

Because of the uncertainty and the complexity of the 
components in the interaction kernel, such as the aver- 
age density of galaxies in the environment n g (t), the rela- 
tive velocity distribution of the aggregation pairs V(t) and 
the interaction cross section E, the exact prescription of 
the interaction kernel is still poorly known. We will thus 
adopt in our simulation a simplified formula for the aggre- 
gation kernel K{m, to , t) with separated time evolution 
term and the cross section term. We give the interaction 
kernel K(m, to , t) as: 



K(m, to , t) oc t k {m 2 ^ ■ 



,2/3 

■ to ) 



1 + a(m 2 / 3 + to 



2/3 



(3) 



where t k is the cosmic time evolution term. The free pa- 
rameter k depends on the specific structures like clusters, 
filaments or sheets; a is a free parameter in our model, 
which describes the relative importance of two kinds of en- 
countering (geometric collisions and focused interactions) 
in the cross section. The aggregation time scale is r, and 
thus the aggregation rate is r , which is proportional to 
the interaction kernel K(m, to , t) as t -1 oc K(m, to , i). 



3. Modelling the ultraluminous infrared phase at 
high redshift 

The binary aggregation dynamics and Monte Carlo sim- 
ulations give us information about the evolution of the 
galaxy mass distribution function N(m,t). What we can 
actually observe and use to constrain the galaxy evolution 
model is: 1) the luminosity functions of galaxies in gen- 
eral or with certain morphologies; 2) the source counts, 
redshift distributions and the background intensity from 
various Hubble Deep Field surveys and the background 
measurements . 

In this section, we will discuss the conversion of the 
mass distribution function N(m, t) by Monte Carlo simu- 
lation of a merger-driven galaxy evolution scenario to the 
observable luminosity function N(L,t). 

In this case, we need to know i) a simple prescription of 
mass-to-light ratio for starburst galaxies and AGNs; and 
ii) the bolometric correction for certain spectral charac- 
teristics, especially at infrared wavelengths. 

A simple discussion of the mass-to-light ratio for the 
starburst galaxies is given by Cavaliere & Menci (1997) 
for the faint blue galaxies who estimate the luminosity of 
starburst galaxies by the gaseous mass of the galaxies and 
the dynamical time scales. The mass to light relation of 
those blue starburst galaxies is given as: 

L ( M \ n 

— =f— J and L*(z) oc f(z, A ,O ) (4) 

where r\ = 4/3; L is the bolometric luminosity and M is 
the mass of the galaxy; L* is the local standard charac- 
teristic luminosity, with the corresponding mass M* (Pee- 
bles 1993). L*(z) oc f(z, Ao,rio) represents a cosmological 
redshift dimming. Considering the large scale structure, 
Cavaliere et al. (1997) gave a prescription of f(z) being a 
function of dimensionality D of the large scale structure 
as L* oc f(z) oc (1 + z)( 3+£> )/ 2 . Studies of the origin of ul- 
traluminous infrared galaxies such as Arp220, NGC1614, 
NGC3256 and IRAS18293-3413, show that the infrared 
luminosity of this kind of galaxies is about a factor of 20 
higher than that of normal starburst galaxies, with a sta- 
tistical relation of Lfi r oc MdT\ (Taniguchi & Ohyama 
1998). The increase in both the dust mass Md and the 
dust temperature Td by starburst merger events could en- 
hance the far infrared luminosities of starburst merging 
galaxies dramatically; we thus call them ultraluminous in- 
frared galaxies (ULIGs). Because of the uncertainties of 
the gas/dust ratio and the complexity of the consequen- 
tial heating of the dust grains, the mass-infrared lumi- 
nosity correlation for the ultraluminous infrared galaxies 
is still unclear. In our calculation, we thus simply adopt 
the power law relation for the mass-to-infrared luminos- 
ity similar to that of starburst galaxies by Cavaliere & 
Menci (1997), with the exponent rj increased by a factor 
of about two to simulate the situation of an enhancement 
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of their infrared luminosity by about a factor of near 20 
for a typical ULIGs of mass ~ 10 12 M Q . 

In our model, we basically assume that starburst and 
AGN activities which are triggered by the merger events 
at high rcdshift are more drastic than those of low redshift, 
simply because the mergers are usually between gas rich 
systems at high redshift and the progenitors of the low 
redshift mergers are already poor in cold gas. In this case, 
we assume in our model a redshift depenent infrared burst 
phase, which means the enhancement of infrared radiation 
from gas rich mergers at high redshift is much higher than 
that of low redshift. We simulate this effect with a power 
law suppression of the infrared burst luminosity simply as 
Li r (z — Az) oc Li r (z)~^ below a transition redshift z ~ 1, 
besides a normal redshift dimming defined in Eq. This 
power law suppression means that the infrared luminous 
galaxies at the bright tail of the luminosity function be- 
come gas poor faster than the less luminous ones. On the 
other hand, choosing redshift z ~ 1 as a redshift transi- 
tion for the gas rich mergers at high redshift and gas poor 
mergers at low redshift is based on the consideration that 
the cosmic time scale of z ~ 1 is about 3 x 10 9 years, 
which is approximately the time scale for the disk evolu- 
tion (Lin & Pringle 1987), it probably indicates a stage 
when galaxies are becoming gas poor. 

While the basic concept is clearly correct, the details 
of these assumptions are a bit crude and arbitrary, espe- 
cially the assumptions for the transition redshift z ~ 1 
and the consequential power law suppression of the in- 
frared burst phase. But it seems that they are the impor- 
tant effects for the fitting of IRAS source counts in our 
simulation, especially the strong evolution at flux range 
10 mJy ~ Uy. Varying any of the model parameters 
could influence more or less the evolution of the luminos- 
ity functions and the source counts finally. While all of 
the model parameters combine to influence the luminos- 
ity functions together, the impact from variation of one 
of the parameters such as a, D or rj, probably still could 
be compensated by varying again other related parame- 
ters. As for the assumptions of transition redshift and the 
differential burst phase, there seems to be not much room 
left for the variation of these particular values. We checked 
the case when the transition redshift and the enhancement 
of the infrared burst phase 77 are both twice and half of 
the values in our model. Although we adjusted other pa- 
rameters in order to approximately fit the observed local 
luminosity functions for the infrared luminous sources, it 
seems that they are far from giving an acceptable fit to 
the IRAS deep surveys (see FigJ3). In this case, what wc 
show here is only one of the most possible, realistic set 
of model parameters which could give the best fitting of 
available observations. Since our model is still too crude, 
wc need probably further information about the evolution 
of infrared bright sources at high redshift, or the redshift 
distributions to high redshift, to give more robust model 
constraints. We will discuss this point again later. 
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Fig. 1. This is our fit to the IRAS 60 /im source counts us- 
ing three major populations (starburst galaxies, Seyferts 
and spirals). The data are from IRAS Point Source Cata- 
log (1985) with galactic latitudes |6| > 50° (PSC), Hacking 
et al. IRAS deep survey (HCH), FSC from deep surveys by 
Moshir et al. (1992) and Saunders (1990). We include here 
also the data by Gregorich et al. (1995) (GNS). It seems 
that the GNS data could overpredict the IR background 
and probably their sample has a biased detection (Bertin 
et al. 1997). The source counts of starburst galaxies and 
Seyferts are from the Monte Carlo simulation where the 
evolution of both activities are triggered by the galaxy- 
galaxy interactions/mergers during the structure forma- 
tion and evolution. The spirals are assumed to have an 
almost constant star formation history since their forma- 
tion. The wiggles in the curve are from the numerical inac- 
curacy of the binning during the Monte Carlo simulation. 
We consider the excess of the data of ~ 10 Jy to be a cos- 
mologically local phenomenon, since we are inside a local 
sheet. 

4. Numerical simulation and discussion 

We use a Monte Carlo inverse-cascading process to simu- 
late the binary aggregation evolution of galaxies which is 
described by Eq. (|l]). In our simulation, we consider the 
initial galaxy mass distribution as a ^-function starting 
from rcdshift z = 15 and of mass M = 2.5 • 1O 1O M0 as an 
approximation. Since binary aggregation dynamics does 
not strongly depend on the initial condition, the memory 
of the initial condition will disappear after the transients 
and the mass distribution is independent of initial details 
with self-similar evolution. The specific discussion of this 
process was presented by Cavaliere et al. (1991, 1992) an- 
alytically and numerically. 

In this case, evolution of the dwarf galaxies with mass 
less than ~ 1O 1O M0 is not relevant to our results. Even 
for those dwarf galaxies of mass ~ 10 10 Mq, we could see 
that their influence is not very important since our re- 
sults only strongly depend on the evolution of infrared- 
luminous galaxies. 
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The large scale structures (LSSs), which have the ad- 
vantage of a larger density and lower velocity dispersions 
compared with the field and the virialized clusters, present 
an ideal environment for galaxy-galaxy interactions to 
take place in. However, they are still quantitatively less 
understood both from the observational and theoretical 
point of view. In the paper about the evolution of faint 
blue galaxies, Cavaliere & Menci (1997) simply discussed 
the relation of the evolution of merger rate which is rep- 
resented by the term t k in Eq. (||) with different envi- 
ronments, such as the homogeneous "field" , the virialized 
clusters or groups and the large scale structures (sheets 
and filaments). They include a free parameter D in the 
expansion of Universe for the Large Scale Structures. 



10mJy<S_{60}<1Jy 



/(*) = 



1 + z 



-i (3+D)/2 



1 + Zi 



(5) 



where Zi n is the redshift when galaxy aggregation becomes 
effective; D is the dimensionality of the large scale struc- 
ture, with D — 2 for shcetlike structures and D = 1 for 
filaments. 

We can transform f{z) to be a function of cosmic time 
t with the conversion of z to t for a flat universe (see 
Peebles 1993) as: 



1 



r 



3 Hofl 1 / 2 



2/3 



(6) 



so, we get fit) oc t k , with k at range of (-4/3, -5/3) corre- 
sponding to D — 1 (filaments) and D = 2 (sheets). 

In our simulation, the best fit structure for the lumi- 
nous infrared galaxies are more like in the filaments with 
D ~ 1 and k ~ —4/3. The free parameter a in equation 
(||) , which describes the relative weight of two kinds of in- 
teractions (geometric collisions or gravitationally focused 
interactions), is close to 0.9. But we know from our simu- 
lation, although varying these parameters could influence 
the final results, we probably still could compensate such 
effects with the adjustment of other relevant model pa- 
rameters. So, what we show here is only one set of the 
most possible values which could give the best fit to the 
present observations. 

In our model, we assumed a redshift dependent in- 
frared burst phase for the merger triggered luminous in- 
frared sources since a transition redshift z ~ 1, which 
means the enhancement of the far infrared luminosity of 
ULIGs from merging galaxies at high redshift is much 
higher than that of low redshift, simply because the high 
redshift mergers are usually between gas rich systems but 
between gas poor systems at low redshift. We found in our 
simulation although the assumption of the "differential 
burst" phase and the consequential power law supression 
is very arbitrary and crude, it is actually very important 
for the fitting of the strong evolution of IRAS 60 fim source 
counts at flux range of 10 mJy ~ 1 Jy. We show the re- 
sults in Fig.0, where the starburst galaxies and Seyferts 
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Fig. 2. This is the redshift distribution of starburst galax- 
ies, Seyferts and spirals at flux density range of 10 mJy ~ 
Uy. The high redshift contribution (z ~ 1.5) is the ul- 
traluminous infrared burst phase from the gas rich merg- 
ers of high redshift. The newest result from the ISO far- 
infrared survey at 175 fim suggests that half of the FIR- 
BACK sources are probably at redshifts greater than 1 
(Dole et al. 1998), which gives further "motivation" for 
considering the model. 



are both assumed to be of the evolution which is driven by 
galaxy-galaxy mergers/interactions, while spirals basically 
keep a constant star formation history since their forma- 
tion. We normalized our Monte Carlo simulation with the 
observed local luminosity function at 60 fim of starburst 
galaxies from Saunders (1990) and Seyferts from Rush et 
al. (1993). Varying this transition redshift by a factor of 
two already has significant effects for the fitting of the 
source counts (see Fig. ||). But, since our model is still too 
crude, we probably only can say what we present here is 
a most probable scenario which gives best fit to the IRAS 
deep surveys and other available observations. 

We consider in our simulation both types of AGNs 
with the assumption that the abundances of Seyfert I and 
Seyfert II at low redshift are approximately equal, which is 
suggested by the extended 12 fim galaxy sample (Rush et 
al. 1993), and the recent Hubble Space Telescope imaging 
survey of nearby AGNs (Malkan et al. 1998). At high red- 
shift (z > 1), we basically assume that the dust shrouded 
phase dominates and has a fraction of about 80% which 
is also suggested by the recent cosmic X-ray background 
(Gilli et al. 1999). 

We choose NGC1068 IR spectra as the standard tem- 
plate for all the obscured AGNs at low redshift, while 
all "Type I" have a SED well represented by the mean 
SED(spectral energy distribution) of Seyfert I. We also 
assume that the early phase of these AGNs show the 
typical spectra such as the dust shrouded F10214+4724, 
and a phase poor in cold gas like the Cloverleaf quasar. 
The template of all these spectra were well modelled by 



6 



Wang et al.: Galaxy evolution and the IR background 



100000 

:p 10000 
m 

r 1000 



m 
cj 

CO 
CO 

1 

"D 



100 



10 



1: - 

2: 

3: 

FSC * 
HCH 
PSC * 


N w% 







1 

0.01 0.1 



10 100 1000 100001 000001 e+06 
SJ60} (mjy) 



Fig. 3. This is the fitting of fRAS 60 fim source counts for 
different transition redshift and the infrared burst phase 
7]. case 1: z = 2, 77 = 4.6; case 2: z = 1, rj = 2.3; case 
3: z = 0.5, ?/ = 1.4. We could see that case 2 is the best 
fit model parameters, meanwhile this set of parameters 
has the best fit also to the local luminosity functions of 
infrared luminous sources. 



Rowan- Robinson (1992), Rowan- Robinson et al. (1993) 
and Granto et al. (1994, 1996, 1997). 

In order to understand the contributing sources for 
the faint slope of the IRAS 60 /im source counts, we 
plot out also the redshift distributions at flux range of 
10 mjy ~ 1 Jy for the three important populations (star- 
burst galaxies, Seyferts and spirals). We see from Fig.|] 
that this faint slope comes from the low redshift starburst 
galaxies which peak at z ~ 0.5, the local spirals with mean 
redshift about 0.1 and the infrared burst phase of high red- 
shift gas rich mergers which peaks at z ~ 1.5. Fig.^ can be 
used to make powerful predictions about the redshifts of 
faint infrared sources. It seems to imply that about two- 
thirds of the faint 60 /im sources should have redshifts 
from a little less than 1 to 2 (and that a fair fraction of 
those will be Seyferts). Recent ISO and NICMOS results 
of the high redshift ultraluminous infrared galaxies further 
increases our motivaton to consider the model (Dole et al. 
1998, Benitez et al. 1998). 

We calculate the background level vl v at 60 /im, it is 
approximately 1.9 nW mT 2 sr^ 1 . We then extrapolate our 
simulation to wavelengths of 25 /im and 100 fim. The back- 
ground intensity are all shown in Fig.^. Clearly, our ex- 
trapolation to 25 [im is too low, compared to other models, 
which suggests that we have too low a mixture of interme- 
diate dust temperatures in our templates for the emission 
spectra. However, our 60 /im background is very close to 
the lower end of the range empirically derived by Malkan 
& Stecker (1998), and so can be considered as a firm lower 
limit. 
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Fig. 4. The extragalactic IR background spectrum pre- 
dicted by various evolution models and from the COBE 
results (Dwek et al. 1998). The DIRBE 2a upper limits at 
1.25 - 60 fim, and detections at 140 and 240 /im with ±2tr 
error bars are represented by solid circles (Hauser et al. 
1998). The FIRAS 125 - 5000 /im detection (Fixsen et al. 

1997) is shown by a light dashed line; the UV-optical lower 
limits of Pozzetti et al. (1997) and the 2.2 /im lower limit 
(Gardner et al. 1997) are represented by solid squares. 
The "X" represent the upper limits on the EBL derived 
from TeV observations of Mrk 501 (Stanev & Franceschini 

1998) , and Mrk 421 (Dwek & Slavin 1994), and the open 
small diamonds represent the upper limits derived from 
the analysis of the fluctuations the DIRBE maps (Kashlin- 
sky et al. 1996). The heavy solid curve represents the EBL 
spectrum calculated from the model with the UV/optical 
indicated SFR and the dust extinction effects (Dwek et al. 
1998). Other models presented in the figure are the Back- 
wards Evolution models of Malkan & Stecker (1998; dot- 
ted line), and of Beichman & Helou (1991; dashed dotted 
line); the Forward Evolution models of Franceschini et al. 
(1997, dashed line), and Guiderdoni et al. (1997, dashed 
triple dotted line); and the Cosmic Chemical Evolution 
model of Fall et al. (1996, thin solid line). The dashed line 
marked S represents the EBL spectrum calculated for a 
dust— free universe with the UVO cosmic star formation 
rate. The lower limits from IRAS counts at 25, 60 and 
100 fim by Hacking & Soifer (1991) are also represented 
here by the solid squares. We plot out the revised lower 
limits at these wavelengths from our model with three big- 
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5. Conclusion 

We described in this paper a Monte Carlo simulation of 
the inverse-cascading process of a merging-driven galaxy 
evolution scenario, where the evolution of infrared lumi- 
nous starburst galaxies and AGNs may be triggered by 
the galaxy-galaxy interactions till even moderate redshift 
(say, z < 1) in the Large Scale Structures. We assume 
in our model a redshift dependent infrared burst phase 
which is based on the concepts that starburst and AGN 
activities triggered by gas rich mergers at high redshift 
are more drastic than those of low redshift, thus the en- 
hancement of the far infrared luminosities in these ULIGs 
from the high redshift merger events are higher than that 
of low redshift. We simulate this effect in our calculation 
by a power law suppression of the infrared burst phase 
since a transition redshift z ~ 1. We adopt the transition 
redshift here at z ~ 1 simply because the cosmic time 
scale of z ~ 1 (3 x 10 9 years) is approximately the disk 
evolution time scale (Lin & Pringle 1987). Varying any of 
the model parameters could influence the evolution of lu- 
minosity functions of infrared luminous sources and thus 
influence the souce counts more or less. But, no matter 
how we adjust the relevant parameters, such as a, D and 
77, in order to get a strong decrease of merger rate with 
cosmic time around redshift 1 ~ 2 for the source count fit- 
ting, it seems that the quick fading and suppression of the 
infrared burst phase at redshift z ~ 1 is a critical effect 
in our simulation to interpret the strong evolution of the 
IRAS 60 /im source counts at flux range 10 mJy ~ 1 Jy. 
The impact from the variation of model parameters such 
as a, D and r), probably could be compensated by vary- 
ing again other related parameters, thus their particular 
values in our model are not the critical points for such 
an evolutionary scenario; However, varying the transition 
redshift and the infrared enhancement r] by a factor of two 
could strongly influence the final results. So, it appears 
that there is not much room for varying any of these pa- 
rameters. Since our model is still too simple and crude, 
we probably need further information about the evolution 
of luminosity functions at high redshift for those infrared 
bright sources in order to give strong model constraints. 

We checked also the redshift distribution of the three 
major infrared sources (starburst galaxies, Seyferts and 
spirals). We see from Fig.|] that the mean redshift of 
starburst galaxies and AGNs which are brighter than 
5*60 ~ lOmJy is around z ~ 0.5 and quickly diminish 
till z ~ 1; a new population of high redshift ultralumi- 
nous infrared burst phase at mean redshift z ~ 1.5 then 
takes over. This seems to be consistent with the present 
groundbased NIR and optical/UV HDF surveys, where 
they failed to have enough detection of starburst galaxies 
near z ~ 1 (Ashby et al. 1992, Koo & Kron 1992). On the 
other hand, recent NICMOS/VLT and FIR/submm sur- 
vey surely found a certain amount of infrared bright galax- 
ies at high redshift, and the newest result of FIRBACK 



(Far Infrared Survey) with ISO shows that more than half 
of the ultraluminous infrared galaxies are at redshift z > 1 
(Lilly et al. 1998, Benitez et al. 1998, Dole et al. 1998). 
This provides further motivation to consider their contri- 
bution to the strong evolution of the IRAS faint source 
counts at flux range of 10 mJy ~ 1 Jy in Figj^. 

Meanwhile the new data of nearby blazars like Mkn501 
by HEGRA team, appears to show an extension of the 
TeV 7 ray spectrum till about 24 TeV, which sets an up- 
per limit for the intergalactic infrared emission history 
during the structure formation (Aharonian et al. 1999). 
We calculated the background level at 60 fxm from a pos- 
sible merger driven starbursts and AGNs scheme and 
a simple constant star formation history for the spiral 
galaxies. The infrared background level at 60 fim is only 
1.9nWm~ 2 sr _1 , which is about half of those estimated 
by some previous papers and consistent with the upper 
limit from the new data of TeV 7 ray spectrum of Mkn501 
(Stecker 1999, Funk et al. 1998). Clearly, this is a strong 
lower limit, because any variation on our model could pro- 
duce a higher background. The forth-coming data on di- 
rect source counts at infrared wavelengths will allow to 
further constrain the evolution of both AGN and star- 
bursts, as well as the absorption at gamma rays near 10 
TeV photon energy. 
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